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Abstract 

Recent analysis shows that, if the 125-126 GeV LHC resonance turns out to be the Standard 
Model Higgs boson, the electroweak vacuum would be a metastable state at 98% C.L. In this 
paper we argue that, during inflation, the electroweak vacuum can actually be very short-lived, 
contrary to the conclusion that follows from the flat spacetime analysis. Namely, in the case of 
a pure Higgs potential the electroweak vacuum decays via the Hawking-Moss transition, which 
has no flat spacetime analogue. As a result, the Higgs vacuum is unstable, unless the rate of 
inflation is low enough: H mi < 10 7 — 10 10 GeV. Models of inflation with such a low rate typically 
predict negligible tensor perturbations in the cosmic microwave background radiation (CMBR). 
This is also true for models in which the perturbations are produced by a curvaton field. We 
also find that if the effective curvature of the Higgs potential at a local maximum (which may 
be induced by inflaton-Higgs interactions) is large enough, then the decay of the electroweak 
vacuum is dominated by the Coleman-de Luccia transition. The electroweak vacuum is also 
short-lived in this case, due to a negative effective self-interaction coupling. Based on our 
analysis of Higgs vacuum stability during inflation, we conclude that the observation of tensor 
perturbations by the Planck satellite would provide strong indirect evidence for new physics 
beyond the Standard Model responsible for stabilisation of the electroweak vacuum. 



1 Introduction 



The discovery of a Higgs-like resonance at 125-126 GeV [T], [2J, triggered much discussion on 
the stability of the electroweak vacuum. In the recent, most accurate analysis [3], it was found 
that, if the Standard Model is valid up to Planckian energies, a Higgs mass < 126 GeV 
implies that the electroweak vacuum is a metastable state at 98% CL. Any discussion of Higgs 
vacuum stability must be considered in a cosmological framework. The classical flat spacetime 
bounce solutions jl], [5], which describe the transition from a false vacuum to a true vacuum 
in the semiclassical approximation, usually have small gravitational corrections [6], [7], and, in 
fact, the bound quoted in [3] is considered in flat spacetime. 

There are strong theoretical arguments as well as indications from astrophysical observations 
in favour of an inflationary period of the early universe, which is well-approximated by a 
cosmological patch of de Sitter spacetime. Inflation not only solves some problems with the 
standard hot Big Bang theory, such as the horizon and flatness problems, but also provides 
a natural mechanism for the generation of large scale density (scalar) and gravitational wave 
(tensor) perturbations, imprinted as anisotropies in the cosmic microwave background radiation 
(CMBR) [8]. False vacuum decay in (quasi-)de Sitter spacetime may differ dramatically from 
that in flat spacetime because there are known de Sitter bounce solutions which have no flat 
spacetime analogue and, moreover, in certain situations they dominate false vacuum decay 
processes. One of these solutions is the Hawking- Moss instanton [9]. It usually provides the 
dominant contribution to false vacuum decay when a scalar potential around its local maximum, 
separating false and true vacua, is flat enough. In this paper we show that this is indeed the 
case for the Higgs potential and find that the electroweak vacuum is unstable during inflation, 
unless the inflationary rate, H ini , is low enough: H inf < 10 7 — 10 10 GeV. We also consider the 
situation in which large curvature of the Higgs potential is induced during inflation by inflaton- 
Higgs interactions. In this case we find that electroweak vacuum decay is dominated by the 
Coleman-de Luccia instanton, and the vacuum turns out to be unacceptably short-lived. 

Interestingly, our analysis shows that important information about electroweak vacuum sta- 
bility can be inferred from a possible observation of tensor fluctuations by the Planck satellite 
[10J. Such an observation would exclude low-scale inflationary models preferred by the require- 
ment of vacuum stability, and would provide a strong indication of new physics responsible for 
stabilisation of the electroweak vacuum. 

The rest of the paper is organised as follows. In the next section we briefly review electroweak 
vacuum stability in flat spacetime. In section 2 we discuss Hawking-Moss and Coleman-de 
Luccia transitions in de Sitter spacetime, followed by a section in which we discuss constraints 
on models of inflation from the requirement of electroweak vacuum stability. 

2 Electroweak vacuum stability in flat spacetime 

The decay rate of the electroweak vacuum is calculated as follows. Since the instability scale 
(which is defined as the scale /i; where the running Higgs quartic coupling A(/z) changes sign from 
positive to negative, that is, A(/ii) = 0) is much higher than the electroweak scale: /i; 3> v^w, 
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Figure 1: The Higgs potential. For large values of the Higgs field h, the electroweak vacuum config- 
uration is regarded as trivial, vew ~ 0. 

one can neglect the Higgs mass parameter and thus approximate the 1-loop Higgs potential as 

m 

Vh = , (1) 

where X c s(h) = A(/x) + {3\ln Here we ignore the numerically less significant quantum 

correction arising from wave-function renormalisation and (3\ is the 1-loop /3-function that 
defines the running of the Higgs quartic coupling, A(/i). Since the potential ([TJ does not 
depend on the renormalisation scale, fi, we find it convenient to evaluate the potential at the 
instability scale, Hence we have: 

Vk = f^-h* , (2) 

where (3\ is the {3\ at \i = //j. 

The stationary equation, = 0, admits two solutions. One is the local minimum, hf v = 0, 
which, in the given approximation (ti EW ps 0), represents the electroweak vacuum configuration. 
Another, h* = /i;e _1//4 , represents a local maximum since (3\ < in the Standard Model and 
thus the curvature of the potential, //*, at h = /i*, is negative: = = jf^/if < 0. 

The potential is sketched in Fig. [T] 

In flat spacetime the transition from the electroweak false vacuum to a deeper (negative 
energy) true vacuum is described by the 0(4)-symmetric Lee- Weinberg bounce [3]. The classical 



h=h, 
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Euclidean action evaluated at this configuration is: 



3|A(/i m )| 

where the renormalisation scale, /i m , (R ~ l/// m is the size of the bounce) is that which 
minimises the action ([3]), i.e. (3\(iJ, m ) = [12J. In the Standard Model /i m ~ 10 17 GeV. Hence, 
very small bounces (i?/ij <C 1) provide the dominant contribution to false vacuum decay and, 
therefore, spacetime curvature effects are expected to be small (this was demonstrated explicitly 
in [?]). Thus the transition probability within the lifetime of the universe, ttj = 1/Hq w 10 32 
eV" 1 , ia3 

p = e -ft» « e 390 -« « 1 , (4) 

since |A(/i m )| ~ 0.01 — 0.02 depending on the top quark mas^J Therefore the electroweak 
vacuum can be regarded as a metastable, but safely long-lived state for m h m 125 — 126 GeV 



3 Electroweak vacuum stability during inflation 

Vacuum stability must obviously be discussed in a cosmological context, where the question of 
initial conditions becomes important. A priori, the magnitude of the Higgs field can assume 
any value in the interval [0, Mp], so the probability of finding it at an initial value less than the 
instability scale, h < ^ ~ 10 10 GeV, is, in fact, extremely low: (/i;/M P ) 4 & 10~ 32 [13]. We 
recognise this as a serious problem, but take, nevertheless, a tolerant attitude. I.e. we assume 
that, as a result of some unknown dynamical mechanism (see, e.g. [14J) or via anthropic 
reasoning, the Higgs field is initially in its false electroweak vacuum state. 

Inflation is a widely accepted paradigm favoured by some astrophysical observations. Ac- 
cepting the inflationary paradigm, we extend our Higgs potential ([I]) by adding terms containing 
the inflaton field <p: 

V(h, 0) = V M {<f>) + V{<f>, h) , where V(0, h) = V K {h) + V K . M {h, 0) (5) 

It is convenient to expand the inflaton potential in a Taylor series around 0i n f, the field 
value just at the start of the observable inflation era, that is, ~ iV = 60 e-folds before the end 
of inflation when the density perturbation on the scale of present Hubble radius was generated: 

v M ((f>) = v inf + - M + \vi'{4> - inf ) 2 + (6) 

1 Actually, e _p gives the probability to find the visible universe in the false vacuum, while (1 — e~ p ) is the 
probability of having the universe in the true vacuum. 

2 The Higgs vacuum stability is very sensitive to the top quark mass. Here we take rat = 173 ± 0.9 GeV 
15j. However, as it has been pointed out in j!6| . uncertainties in top quark mass measurements may actually 
be larger. 
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where V inf = V(0 inf ), K = dV/d^^, V? = d 2 V/d<p : 

0, 4>) = without loss of generality. The potential ^ must be flat enough around 
0), i.e., 

V 



, etc. We assume that V(h 

= 0inf , h 



e 



My, 



V inf 



< 1 , - 1 < 77 



V" 

M 2 -^- < 1 

Vinf 



(7) 



so that the energy density in the universe is dominated by the constant term: V in f ~ Vi n f. This 
leads to a de Sitter-like exponential expansion of the universe characterised by the expansion 
rate H m { = y^j^r, where Mp w 2.44- 10 18 GeV is the reduced Planck mass. Models of inflation 

are broadly classified as large-field ( 2e > r\ > 0), small-field (e > 0, r\ < 0) and hybrid (r/ > 2e). 

To discuss instanton mediated transitions between false and true vacua we consider the 
Euclidean analogue of cosmological spacetime: 



ds 2 = d x 2 + p 2 (x)dn 



2 

3 > 



where dQ$ is the length element on a 3-sphere, p(x) is the Euclidean scale factor and x 2 — t 2 +r 2 . 
The relevant Euclidean action for homogeneous fields is: 
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where A = dA/dx and A = d 2 A/dx 2 - The resulting equations of motion read: 

i + 3 1 - 



(9) 



P 



1 + 



3M2' 





«9V 


p 




h= 
p 


av 


2 






(10) 

(11) 

(12) 



We take <fi ~ 0^ = const, as an approximate solution to Eq. (10) during inflation. Assuming 
that the inflaton provides the dominant contribution to the energy of the system, V ~ Vi n f 3> h 2 , 
it is easy to verify that Eq. (12) admits the solution p = sin (i? in fx) (% G [0, 7r/if in f]). 
Hence, in this fixed background approximation we are left with the following equation for h: 



h + 3Finf COt (#infX) h 



dV(h, 



>infj 



dh 



(13) 



Since cot (H ini x) diverges at x — {0, vr/if inf }, we supplement equation (13) with the boundary 
conditions: h(0) = h{ix / H ini ) = 0. 
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3.1 Vacuum decay via the Hawking-Moss instanton 

Eq. (13) admits the trivial solution, h = 0, which describes the electroweak false vacuum state 
during inflation. Noting that dv ^ uh ^ = Q, we find another simple solution: h = h*. 



ah h=h, 

The transition between the electroweak vacuum state (0; n f , h = 0) and the state (0; n f , 
corresponding to the top of the barrier in the /i-direction (see Fig. [l| during inflation is 
determined by the difference between the corresponding Euclidean actions: 

Bum = S[(f> = inf , h = h*} - S[<f> = inf , h = 0} = (V(0 in f , K) - V inf ) . (14) 

This is the Hawking-Moss transition, which typically dominates in false vacuum decay processes 
when the curvature of the potential near it's local maximum is smaller than the expansion 
rate during inflation [9]. Physically, the Hawking-Moss transition can be understood as a 
thermally-driven transition where the false vacuum state, h = 0, is excited to unstable state, 
h = within the Hubble radius [17]. The de Sitter spacetime plays the role of a thermal 
reservoir with associated Hawking- Gibbons temperature T^s = -£/mf/27r [T7j. Once h = 
there is unit probability that the configuration will subsequently roll down towards the negative 
energy true vacuum state. This, in fact, causes a disastrous collapse of the universe into 
a black hole singularity ('Big Crunch') rather than relaxation into the true vacuum. The 
transition probability from the electroweak false vacuum (0i n f, 0) to a state (0mf , K) during 
'visible' inflation, N = r in{ H in{ w 60, is given by: 



v ; i 1Q7 87T 2 V H (K) + ^nf-H(0infA, , 

p = N e HM w exp <( 13.7 \ (15) 

3 -"inf 



3.2 Vacuum decay via the Coleman-de Luccia instanton 

Besides the unknown inflation rate, H in f, the above probability also depends on unknown 
inflaton-Higgs interactions described by Vinf-H- It is reasonable to require that these inter- 
actions do not spoil the flatness of the inflaton potential, Vinf) m order to satisfy the slow-roll 
conditions <Q, as well as the flatness of the Higgs potential, Vh, in order to keep the Higgs mass 
light. However, if we introduce a large bare mass for the Higgs boson and accept an unnatural 
fine-tuning of parameters, the effective curvature of the Higgs potential near the local maximum 
can be large as a result of inflaton-Higgs interactions. In that case the Hawking-Moss transition 
is suppressed and false vacuum decay via the Coleman-de Luccia instanton becomes dominant. 
To make our discussion quantitative let us consider the following inflaton-Higgs interaction 
potential: Vi n f_H = ^f(4>)h 2 , where /(</>) is an arbitrary positive function of the inflaton field 
0. In the fixed background approximation (0 = const.) this interaction potential turns into 
an effective mass term for the Higgs field, with m 2 cS = }\4>- xn {). Then the total effective Higgs 
potential during inflation can be written as: 



V{(f)- mi , h) = ^m 2 cS h 2 



1 / /i x 2 



2 V h 



(16) 
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where h* is a local maximum, 



m 



cff 



A 



(17) 



for m^jj > 0, and recall that for large values of h we have A < (the potential is unbounded 
from below, see Fig. [I]) so here we approximate it as a constant. 

We are interested in the case of large effective curvature of the potential at the local max- 



imum, that is, m^ ff ^> \Hf ni . In this limit Eq. (13) takes the form of the corresponding flat 
spacetime equation: 



dVj^h) 
dh 



h" + -ti 
x 

where A' = dA/dx, A" = d 2 A/dx 2 and x = m e sX- Decay of the electroweak false vacuum 
proceeds through the nucleation of a bubble, with the field h(x) approaching hn > h* near the 
centre and ht < h* far outside the bubble (see Fig. [TJ. Note that, contrary to the similar process 
of the bubble nucleation in flat spacetime, h L /h R need not necessarily be the false/true vacuum 



in de Sitter spacetime [17J. For the interior bubble solution we approximate the potential (16) 
by neglecting the quartic coupling, while for the exterior solution we neglect the quadratic 
coupling instead. These two solutions must be continuous at x*\ 



(19) 



The following one-parameter class of bounce solutions can be obtained analytically in the 
above approximation: 



h(x) 



x(Ji(ix t )+iY\(—ix,)) 



(Ji(ix) + iY\{—ix)) , < x < oo 



2^2K fh R 

X* = — : I 1 



h R \h 



1/2 



Then we can compute the exponential B-factor: 



CdL 



where 



I — I x 3 dx 
'o 



h 2 {x) 



i h 2 (x) 
2hl 



(20) 
(21) 

(22) 
(23) 



4 Constraints on inflation from vacuum stability 

First we consider the case of negligible inflaton-Higgs interactions, i.e we set Mnf-H = in (J5|. 
Then Vh(^*) is given by (J2l), where h* = /ije -1 / 4 . The condition of vacuum (meta) stability, 
p < 1, implies the following constraint on the rate of inflationary expansion: 

H ini < 4.6 ■ 10 7 (2.7 ■ 10 10 ) GeV , (24) 
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for rrih = 126 GeV and m t = 174(172) GeV. For mh = 125 GeV the limits are an order of 
magnitude lower. Assuming that quantum fluctuations of the inflaton field are entirely respon- 
sible for the production of density perturbations that are imprinted in the CMBR temperature 
anisotropy, we can compute the standard relations between various observables and the slow-roll 
parameters ([8]). In the linear approximation they are: 

As = (8^r V2 ff , (25) 
n s = 1 - 6e + 2r] , (26) 
r = 16e , (27) 

where A s is the amplitude of scalar perturbations, n s is the spectral index and r is the ratio 
of squares of scalar and tensor perturbation amplitudes. Current experimental values for the 
above observables are [19J: 

A c * p ~ 5 • 10" 5 , (28) 
n fP = 0.9538 ± 0.0081 , (29) 
r exp < 0.11 , (30) 



As follows from (25), in order to reproduce the observed amplitude of density perturbations 
(28) with the inflation scale (24), one must have e « 10~ 16 (10~ 10 ). The observed negative tilt 



(n s < 1) in the spectral index (29) then implies rj < and, therefore, only small-field models are 
viable. Note that the tensor-to-scalar ratio is also tiny in this case: r « 10 _15 (10~ 9 ), well below 
the sensitivity (~ 10~ 2 ) of the Planck satellite experiment JTDJ to CMBR B-modes. Therefore an 
observation of tensor perturbations by Planck would rule out small-field inflationary inflationary 
models, leaving the problem of vacuum stability unresolved within the Standard Model. 

As discussed in the previous section, sizeable inflaton-Higgs interactions, VJ n f_H, may in- 
crease the effective curvature of the scalar potential near the local maximum, h = /i*, during 



inflation and thus reduce the Hawking-Moss transition probability (15). This possibility was 
recently discussed in [H] as a mechanism to ensure vacuum stability during inflation. Following 
[14] consider, for example, V| n f_H = ^h 2 (p 2 (with a > 0). Then the effective mass, m e s, during 
inflation, is given by = a0f nf and = (— a/A) 1 ^ 2 in f > /iie" 1 / 4 (\(h*) < 0), so that the 



transition probability (15) is suppressed by the large negative factor in the exponent. This can 



be achieved in large-field (0; n f > Mp) chaotic inflation models providing 



a > 1.4vffl ( • (31) 



It is interesting that current data [18] . [19] disfavors all large-field models except the simplest 
one with quadratic potential^ VJ n f = |m|0 2 , where the inflaton mass is ~ 10~ 5 M P . Taking 



3 The model with a quadratic potential can also be viewed as an approximation of the model of natural 
inflation VM. 
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|A| ~ 0.01, we obtain a > 6.0 • 10~ 12 within this model. It was estimated in [H] that, to 
have a sufficiently flat quantum-corrected inflaton potential, one needs to assume a < 10 -3 . 
However, quantum corrections from inflaton-Higgs interactions also give a contribution to the 
Higgs potential, i.e. the Higgs mass receives a finite radiative correction, which in one-loop 
approximation reads: 

Sm l « ^~2 m l ■ ( 32 ) 

By demanding that this correction not exceed the observed Higgs mass: 5rrih < rrih, we obtain 
a < 647r 2 (mft,/m^) 2 ~ 1.7- 10~ 20 , which is incompatible with the above vacuum (meta) stability 
condition ( |3lj ). Note that a very small a parameter is technically natural, because in the limit 
a — > the Higgs and inflaton sectors decouple (see, for example, the discussion in [2T]). 

However, if we abandon the naturalness condition, we can always introduce a bare Higgs 
mass parameter and tune it with radiative corrections to obtain the observed Higgs mass 
without further restriction on the strength of the inflaton-Higgs interaction. In this case the 
contribution to false vacuum decay is dominated by the Coleman-de Luccia transition described 



by Eq. (22). Numerically, we find that the integral (23) is dominated by a negative contribution 
from the interior solution, i.e., 7 < 0. Since the self-interaction coupling is also negative, the 
bubble nucleation rate is p ~ e~ BcdL 3> 1, and thus the electroweak vacuum is again unstable. 

Finally, there is also the possibility that CMBR perturbations are produced by quantum 
fluctuations of a field other than the inflaton field, as in the curvaton scenario [22]. Then the 



observables (25)-(27) are modified to 



^• = ^ + ! W' (33) 
< = 1 - 2£ -rr^)' (34) 
r ' = rrfrn • < 35 > 

where g 2 (X) = ^ for a 'large' initial curvaton field value, a*, with X = a^/Mp ^> 1, and 
g 2 {X) = J? for a 'small' initial curvaton field: X <g 1. With the above equations we see that, 
if the curvaton is a major source of density perturbations then g w 10 5 — 10 8 is required for 
low scale inflation models. This again implies vanishing tensor perturbations: r' w 0. Thus we 
come to the universal conclusion based on the above vacuum stability analysis: an observation 
of tensor perturbations using the upcoming data collected by the Planck satellite would provide 
a strong hint of new physics responsible for stabilisation of the electroweak vacuum. 

5 Conclusion 

We have considered the question of electroweak vacuum stability in an inflationary universe. 
We found that electroweak vacuum decay is governed by the Hawking-Moss transition and it is 
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unacceptably fast, unless the rate of inflation is H in { < 10 7 — 10 10 . We also considered inflaton- 
Higgs interactions which may induce large curvature in the effective Higgs potential during 
inflation, in which case vacuum decay is dominated by the Coleman-de Luccia transition. We 
have demonstrated that the electroweak vacuum is also unstable in this case, essentially due to 
the negative effective Higgs self-interaction coupling. 

These observations have an interesting implication for CMBR fluctuations. Models of infla- 
tion with a low inflationary rate typically predict very small amplitude tensor fluctuations, with 
the same conclusion following from curvaton models. Thus, the upcoming data from the Planck 
satellite may provide important information concerning stability of the electroweak vacuum. 
Namely, the observation of tensor fluctuations would provide strong indirect evidence in favour 
of new physics beyond the Standard Model, responsible for stabilisation of the electroweak 
vacuum. 

Acknowledgments. This work was partially supported by the Australian Research Council. 

References 

G. Aad et al. [ATLAS Collaboration], [arXiv:1207.7214| [hep-ex]]. 



[2] 
[3] 

[4 

[5] 
[6] 
[7 



[9 
[10 

[11 
[12 



S. Chatrchyan et al. [CMS Collaboration], arXiv: 1207.7235 [hep-ex]]. 



G. Degrassi, S. Di Vita, J. Elias-Miro, J. R. Espinosa, G. F. Giudice, G. Isidori and 
A. Strumia, JHEP 1208, 098 (2012) jarXiv: 1205. 6497] [hep-ph]]. 

S. R. Coleman, Phys. Rev. D 15, 2929 (1977) [Erratum-ibid. D 16, 1248 (1977)]; 
C. G. Callan, Jr. and S. R. Coleman, Phys. Rev. D 16, 1762 (1977). 

K. -M. Lee and E. J. Weinberg, Nucl. Phys. B 267, 181 (1986). 

S. R. Coleman and F. De Luccia, Phys. Rev. D 21, 3305 (1980). 

G. Isidori, V. S . Rychkov, A. Strumia and N. Tetradis, Phys. Rev. D 77, 025034 (2008) 
|arXiv:0712.0242 [hep-ph]]; S. Kanno, M. Sasaki and J. Soda, Class. Quant. Grav. 29, 
075010 (2012) [arXiv: 1201.2272| [hep-th]]. 

V. F. Mukhanov, H. A. Feldman and R. H. Brandenberger, Phys. Rept. 215, 203 (1992). 
S. W. Hawking and I. G. Moss, Phys. Lett. B 110, 35 (1982). 



[Planck Collaboration], astro-ph/0604069 



G. Altarelli and G. Isidori, Phys. Lett. B 337, 141 (1994). 

G. Isidori, G. Ridolfi and A. Strumia, Nucl. Phys. B 609, 387 (2001) [hep-ph/ 0104016] . 



9 



M. P. Hertzberg, |arXiv:1210.36~2~4| [hep-ph]. 



O. Lebedev and A. Westphal, arXiv:1210.6987 [hep-ph]. 



[Tevatron Electroweak Working Group and CDF and DO Collaborations], arXiv:1107.5255 
[hep-ex] . 



S. Alekhin, A. Djouadi and S. Moch, Phys. Lett. B 716, 214 (2012) |arXiv:1207.0980 
[hep-ph]]. 



A. R. Brown and E. J. Weinberg, Phys. Rev. D 76, 064003 (2007) |arXiv:0706. 1573 
[hep-th]]. 

E. Komatsu e t al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) 
[arXiv:1001.4538| [astro-ph.CQ]]. 

K. T. Story, C. L. Reichardt, Z. Hou, R. Keisler, K. A. Aird, B. A. Benson, L. E. Bleem 
and J. E. Carlstrom et ai, arXiv:1210.7231 [astro-ph.CO]. 

K. Freese, J. A. Frieman and A. V. Olinto, Phys. Rev. Lett. 65, 3233 (1990). 

R. Foot, A. Kobakhidze , K. L. McDonald and R. R. Volkas, Phys. Rev. D 77, 035006 
(2008) [arXiv:0709.2750| [hep-ph]]. 

[22] K. Enqvist and M. S. Sloth, Nucl. Phys. B 626, 395 (2002) |hep-ph/0109214] ; 



D. H. Lyth and D. Wands, Phys. Lett. B 524, 5 (2002) [hep-ph/01 10 002]; 

T. Mor oi and T. Takaha shi, Phys. Lett. B 522, 215 (2001) [Erratum-ibid. B 539, 303 
(2002)] |hep-ph/0110096] . 

[23] T. Moroi and T. Takahashi, Phys. Rev. D 72, 023505 (2005) |astro-ph/050"5339] . 



10 



